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Abstract 

The energy levels of an impurity center in a deep quantum well of 
width L and depth g are studied analytically . Renormalised pertur- 
bative series are constructed in the regions gL 2 << 1 and gL 2 » 1. 
Maximal binding energy and wave function deformation to a quasi- 
twodimensional function are found to occur at a certain L c satisfying 
yJ~g~L c ~ 1. Similar results may be obtained for the impurity in a 
quantum wire, in a dot or in a multiwell structure. 



1 Introduction 



The energy levels of donor and acceptor centers in quantum well, wires and 
dots have been widely studied recently in the literature [1]. 

In the present paper, we show that very simple analytic estimates to the 
energy levels may be obtained from perturbation theory. The perturbative 
series are renormalised to account for large orders corrections or, said in a 
different way, the renormalised series are forced to interpolate between the 
different expansions obtained in different regimes. 

In a deep quantum well of deepth g and width L (in effective atomic 
units), we construct renormalised series for gL 2 << 1, and gL 2 >> 1. The 
starting points are the perturbative expansions obtained in the following 
limits: 

i) gL 2 << 1, gL << 1. Shallow level, three dimensional coulombic wave 

function. 

ii) gL 2 << 1, gL >> 1. Shallow level, quasi two-dimensional wave function 

out of the well. 

hi) gL 2 >> 1, L « 1. Deep level, quasi two-dimensional wave function 
inside the well. 

iv) gL 2 » 1, L >> 1. Deep level, three dimensional coulombic wave 
function. 

At high values of g, the curve E b (L), i.e. the binding energy as a function 
of L, contains the four regimes i) - iv). Indeed, let us consider, for example, 
g — 60. Then, when L < 1/60 the conditions for i) are fulfilled. In the 
interval 1/60 < L < 1/a/60, we may use an approximation like ii). When 
l/y/60 < L < 1, the regime iii) holds and, finally, at L > 1, the wave function 
is certainly as indicated in iv). 

The transition from ii) to iii) is found to ocur at ^fg~L ~ 1, and corre- 
sponds to a maximum of the binding energy. 

2 Impurity in a quantum well 

We start from the effective mass hamiltonian in effective atomic units 
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H = --A--+ge(z), (1) 

where the energy unit is me 4 / (h 2 k 2 ) , m is the electron effective mass, k- the 
relative dielectric constant, g = Vq / [me 4 / (h 2 k 2 )] is the well depth, and 9 is 
the step function 

— J 0' ^1 < ^ < ^2, ^2 — ^i = ^, / 9 x 

Ul2j _ \ 1, outside. [Z) 

Typically, g > 20. From the eigenvalue of H, E, the binding energy is 
defined as 



E b = E w — E, (3) 

where E w is the threshold for the continuous spectrum in the well, i.e. the 
lowest eigenvalue of the one-dimensional hamiltonian 

H w = -\^ 2 +gQ{z). (4) 

It is found from E w = k 2 /2, where k is the smallest solution of 



gL 2 

The asymptotic expressions for E w are the following 



E w = 5 {l-^L 2 + ^ 2 L 4 + ...}, gL 2 «l (6) 

Ew = h{\-^ + g^-( 2 + h)(-gh) + -}« 

gL 2 » 1. (7) 

We shall obtain analytic estimates to E to compute the binding energy 
as a function of L. The parameter g is assumed to be large. 

The exact quantum number of the present problem is the z-projection of 
angular momentum, which will be called J. Without loss of generality, we 
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may consider only positive values of J. To avoid using degenerate pertur- 
bation theory, we will study the first level with a given J. when L — > or 
L — > oo, we recover the unbounded (3D) coulomb problem, and the states 
we study have angular momentum equal to J. 



3 The renormalised series for gL 2 « 1 

When gL 2 << 1, the levels are shallow, i.e. located near the top of the 
barrier, and the wave function is mainly outside the well. The characteristic 
confinement distance along the z-direction is l z , where Z" 1 = \j2(g — E w ) = 
gL + . . .. When gL « 1, the wave function is basically the 3D coulombic 
wave function, (f> 3D . On the other hand, when gL » 1 (keeping gL 2 « 1, 
i.e. l z » L) the confinement length is much less than the Bohr radius 
(which in our units is one), and the wave function is written approximately 
as exp— z/l z <p2D, where 02D is the 2D coulombic wave function. 
Accordingly, when gL << 1 we write 

H = H 3D + g + V u (8) 

where V\ = —gQ*(z), and the complementary step function is defined from 
Q(z) + Q*(z) = 1. V\ will be treated as perturbation. From first order 
perturbation theory, we obtain 



1 gL 
9 ~ 2(J + 1) 2 ~ 2 J + 1 J!(J+ 1) J+1 

g L 3 fa 2 + ^ 2 + q 2 ) 

+ j\(j + iy^ + ° {Lh (9) 

where we have written, z\ = q\L, Z2 = &L, ^2 — <?i = 1, and without loss of 
generality, we will assume that <ri < 0, ^ > 0. The leading contribution of 
second order perturbation theory will be —0.002 g 2 L 2 for the ground state 
(J = 0), and practically zero for J > 0. 

Thus, groupping (6) and (9), we obtain for the binding energy at gL « 1, 

E b = a + ai (gL) + a 2 (gL) 2 + a 3 (gL) 3 + . . . , (10) 
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where a = 1/(2(J + l) 2 ), a x = 1/(2 J+1 J!(J + 1) J+1 ), a 2 = -0.498 when 
J = and a 2 = -1/2 when J > 1, a 3 = -{4 + q ± q 2 + q 2 )/{g 2 J!( J+ 1) J+4 2 J ), 
etc. 

Now, let us consider the opposite limit, gL » 1. The hamiltonian will 
be written as 



H = H 2D + H W + V 2 , (11) 

where V 2 = 1/p — 1/ V p 2 + z 2 will be considered as perturbation, p is the 
polar coordinate in the plane. In first order perturbation theory, we obtain 

E = ~2(J+l/2) 2+Ew ~ W + (12) 

where the first correction is &i = —8 for the g.s. (J = 0), and b 2 = — 2(2J — 
2)!/(2J+ 1)!/(J+ 1/2) 3 for the excited states (J > 1). Consequently, for the 
binding energy, we obtain 

E b = J^- k +b + ..., (13) 

where 6o — (J + l/2)~ 2 /2. Notice that Eq. (13) may be applied to the g.s. 
when gL > 8. If we recall that gL 2 << 1, then g is forced to be greater than 
64. 

Note that the first three terms of (10) and the first two terms of (13) 
show that Ef, is, in certain approximation, "universal", in the sense that 
it depends only on the variable gL. If we consider corrections such as the 
term a 3 (gL) 3 of (10), this universality is lost. Notice also that the impurity 
position appears for the first time precisely in a 3 . 

Once we have the correct behaviour of E b for large values of gL, we may 
construct a "renormalised" series from (10) such that the large orders of 
this series will account for the correct asymptotics at gL >> 1. We follow 
the idea of paper [2] in which the method was applied to the two-electron 
problem in a quantum dot. 

For simplicity, we consider the g.s. of the centered impurity and include 
the first two terms of (10). We write f3 = gL/(a + gL). When gL — > 0, 
/3 — > 0, whereas when gL — > oo, /3 — > 1. a is a free parameter which will 
be used to fit the numerical results. It gives an idea of where the transition 
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from "small" gL to "large" gL takes place. The renormalised series is looked 
for as a series in f3 



E b = c + Cl p + c 2 (3 2 + c 3/ 3 3 + .... 



(14) 



The first coeficients c k are defined in such a way that when (3 — > the 
first terms of (10) are reproduced. That is: Co = 1/2, c\ = a/2. On the other 
hand, c 2 and c 3 are chosen to reproduce the correct behaviour at (3 — > 1. We 
note that, as (3 — > 1 

£ b = (c + ci + c 2 + c 3 ) - (ci + 2c 2 + 3c 3 )(l -/?) + ..., (15) 

that is, c 3 and c 4 are required to satisfy: l/2 + o;/2 + c 2 + c 3 = 2, a/2 + 2c 2 + 
3c 3 — 8/ a. 

We show in Fig. 1 variational computations (points) for the g.s. of the 
centered impurity at g — 60, 100, 140 and 180. At each g, values of L 
observing L < 1/ ^fg were included. A trial function with two nonlinear 
parameters was used in the computations [3] . The results show a very small 
dispersion of the points around a fixed curve, in accordance with the predicted 
universality. The curve is well fitted by the series (14) with a ^ A (the solid 
line) . 

4 The renormalised series for gL 2 » 1 

Now, we consider the situation in which the level is deep inside the well, 
gL 2 >> 1. At large values of g, we have again two limiting situations: L « 
1 = a B , and L » 1. The corresponding wave functions are approximately 
(p 2 D sin [71- (2 -zi)/L], and <p 3D . 

When L » 1, we take V 3 = gQ(z) as a perturbation to H 3D . Assuming 
that Z\ and z 2 are both finite, we obtain 



where q is the minimum between and q 2 - The binding energy is thus 




(16) 




(17) 
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where b' Q = ( J + l)~ 2 /2, b' 2 = n 2 /2, b' 3 = -n 2 ^2/~g, etc. 

On the other hand, when L « 1 we write an expression like (11), and 
consider Vi as a perturbation. To first order, we get 

^-2ijTW + E --" ltl + '- (18) 

where the first correction is a[ = — 16[l/2+q+^ — (1— cos 27K 1 )/(27r 2 )] for the 
g.s. (J = 0), and a' 2 = -4(2J-2)!(q 2 + q + l/3-l/2/7r 2 )/(2J+l)!/(J+l/2) 3 
for the excited states (J > 1). That is, 

E b = a' + a' k L k + . . . , (19) 

where a' = 1/2/ ( J + 1/2) 2 . 

Notice that the series (17) and (19) suggest, again, a universal behaviour 
of E b in the leading approximation. That is, E b will not depend on g when 
gL 2 » I. 

For simplicity, we consider again the g.s. of the centered impurity and 
include terms up to 1/L 2 in (17). To construct the renormalised series for 
gL 2 » 1, we write f3' = L/(a' + L), 



E b = c' + 4(1 - f3') + c' 2 (l - (3') 2 + 4(1 - + 4(1 - pr + ■■■■ (20) 

The coefficients c' , c[ and c' 2 are obtained from (17), i.e. c' = 1/2, c[ = 0, 

4 = tt 2 /2/q; /2 . The coefficients 4 an d 4 are required to satisfy the linear 
equations 1/2 + ir 2 /2/a' 2 + 4 + 4 = 2, vr 2 /a' 2 + 34 + 44 = a' {A - 16/vr 2 ). 

We show in Fig. 2 how close the behaviour of E b (L) is to the universal 
behaviour when g = 60, 100, 140 and 180. Variational computations cor- 
responding to values of L for which L > 2.6/ ^fg are represented as points. 
The dispersion of the points is very small. The universal curve is fitted well 
by the series (20) when the parameter a' is near to 1. The result of the fit is 
presented as a solid line. 

5 The maximum of the curve Eb vs L 

The conclusion of the previous analysis is the following. As L is increased 
from to 1/ ^fg the wave function changes from a tridimensional coulombic 
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wave function to a quasi-twodimensional function out of the well. The bind- 
ing energy is a universal function of gL. On the other hand, as L is decreased 
from infinity to 2.6/ \fg the wave function undergoes a change from the tridi- 
mensional coulombic wave function to a quasi-twodimensional function inside 
the well. The binding energy turns out to be a universal function of L. As a 
result of compressing the wave function, the binding energy increases. There 
is a critical value, L c , at which the wave function is maximally compressed 
and Eb is maximal. As indicated, 

1 2 6 

4= < l < < -T- ( 21 ) 

We may obtain a rough estimate to L c as the point at which the two 
series (13) and (19) coincide. For the g.s. of the centered impurity, we get 

L c » — . (22) 

Numerical computations at g > 60 show that the product ^J~g~L c is almost 
constant, taking a value near 1.5 (Fig. 3). 

We note that the dependence L c ~ t/y/g comes also from very simple 
reasonings based on the indeterminacy relations [4]. 



6 Concluding Remarks 

We have shown that, at large values of g, the parameters gL and gL 2 identify 
different regimes in the behaviour of the binding energy and the wave function 
of the impurity problem in a quantum well. In particular, the binding energy 
was shown to be a "universal" function in each of the regions gL 2 « 1 and 
gL 2 » 1. 

Our description is qualitatively valid in problems where g = 20 — 60. As 
g is still decreased, regimes ii) and iii) become less applicable, and may be 
absent at all . In particular, at g ~ 1, we expect a transition directly from i) 
to iv) as L goes from L<<ltoL>>l. 

We would like to stress that there are similarities between the results of 
our section 4 and paper [5], in which the g — > oo limit is studied. In that 
paper, a decomposition like (11) is choosen in the entire interval < L < oo, 
and perturbation theory is applied. The failure of this decomposition at 
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L > 1 is corrected by introducing a scaling parameter and asking for the 
virial theorem to hold. Our approach is more qualitative but, at the same 
time, more exact, stressing what the actual wave function really is at any L. 

A qualitative analysis of wave functions and energy curves for the im- 
purity problem in dots, wires and multiwell structures is as simple as the 
analysis presented above. It may be a useful complement to the existing (or 
in progress) sophisticated numerical calculations. 
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Figure Captions 



Fig. 1. Twice the binding energy of the centered impurity as a function of 
gL. Points are variational calculations satisfying L < 1/^fg. 

Fig. 2. Twice the binding energy as a function of L when L > 2.6/ v /g. 

Fig. 3. The dependence of y/g~L c on g. 
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